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Abstract.In the paper discusses the interaction between two charged 
balls in equilibrium state. It is shown that, depending of the sizes, 
charges and distance, the balls can move in the same or opposite direc¬ 
tion. They can repulse and attract. It is proved, that one of the balls 
may vibrate. 


1. Introduction 

S. Poisson was the first who investigate the problem of determining 
the electrostatic force between two charged balls. Later J. Maxwell [1] 
proved, that the force of interaction between two balls with charges Q 
and g, in eqnilibrinm state, the force different from the force of inter¬ 
action between point charges Q and q located at the centers of those 
balls. The bi-spherical coordinates [2] are used to derive exact formu¬ 
las for the force of interaction between the balls. But the results are 
complicated and inconvenient. Some articles 0, a. 0.0 are dedi¬ 
cated to the task of finding simpler formulas. In this paper we use the 
Kelvin transform and derive formulas for the equilibrium distribution 
of charges on each sphere. The obtained formulas have a relatively 
simple form and permit us to determine the distribution of charges on 
the balls’ surfaces. We also observe certain new effects: it turns out 
that two balls with charges of the same sign may move in the same 
direction and at the same time repulse and attract. In some conditions 
one of the balls can vibrate. 

2. Location of’’ghost - charges” 

We will consider conductors of the form 

E = B{xo,R) U B{yo,r), 
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where 


B{a,R) = {x\ ||a; — a|| </?}. 

Assume that the ball R) carry the charge Q and the ball 
carry the charge q. Denote by 


d = ||fo - foil > R + r. 


the distance between the centers of the balls. 
Let us define Xi by the formula 


Xi= Xo + (fo - ^o) 




Ilfo -^olP 


The point xi is symmetric to fo with respect to the sphere dB{xo, R). 
Similarly, the point 

J.2 

fi = fo + (^0 - Vo)-^ -^ 

l|a^o -yolr 

is symmetric to Xq with respect to the sphere dB{yo,r). 

The points Xn and are the location of ’’ghost - charges”. Note, 
that 

XneB{xo,R), yneB{yo,r), n = 0,1,2,.... 

At the points Xn, iJn let us put charges Qn, dn- The ’’ghost - charges” 
itself do not appear in equilibrium distribution, but out of conductor 
E, generate the same potential function as equilibrium distribution. 

Since the points fm-i, Tm, n = 1, 2,... are symmetric with respect 
to the sphere dB{xo, R) so 


1 


\x 


I/m—11 


||fo-fm-l|| 1 


R 



X G dB{xo, R)- 


Let us introduce two sequences of nonnegative numbers u{n) and 
v{n) by the formulas 


u{n) = \\xn - xo\\, v{n) = WiJn - yo\\, n = 0,l,2,... 
We have u(0) = u(0) =0 and 


u{n) = 


-, v{n) = 


R^ 


d — v{n — 1)’’ d — u{n — 1) 

Let us note that the following inequalities hold: 

r 


, n = 1,2,... 


R^ 


<r, n = l,2,..., 
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3. The ’’ghost - charges” 


The ’’ghost - charges” are 


^ ^ Qn^ (3^ ^n) “I” ^ ^ Qn^ ( 


X-Vn 


For an arbitrary numbers Qq, go dehne the sequences Qn, Qn'- n = 
1,2,... by the following recursive formulas 

_ v{n) 2 

7~j yn—1? yn 'Vn—1? ±, z,, . . . . 

K r 

Let for some M < oo the inequalities IQol, |go| < dL hold. Then 


\q2m\, 

These inequalities we can write in the form 


IQnP, |gnP < 


, m = 1,2,3,.... 


, n = 0,1,2,.... 


Since Rr < {d — R){d — r) so the sequences Qn, dn tend to zero as 
geomteric progression. 


4. The Potential Function of the equilibrium distribution 


Let us denote 


U{x) = J2 


s II* - 


-, X ^ E. 


The total charge Q of the ball B{xo, R) equals: 

^ o q , “(2)*(1) , 1 /'“(l) , ■41H3H2) 

o = = 1 + —7^ + ... -*Hr + —Sv— 


Similarly, the total charge g of the ball B{yQ,r) equals: 


It is well-known that the potential function of the equilibrium distri¬ 
bution U{x) uniquely determined by the following properties 

1 . 

AU{x) = 0 X ^ E] 


U{x) = const, X G B{xo,R); 
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2 . 



3 . 


U{x) = const, X G B{yQ,r)] 


lim U{x) = 0. 


It is obvious that the function (1) satisfies conditions 1, and 4. 
We have 


U{x) = j-^ + qo 
\\x - Xoll 




1 _ ||fn+l - ^oll 

\x-yn\\ R\\x-Xn+l\ 


\\Xl - XqW 

R\\x — xi\ 


= —, xedB{xQ,R). 

K 


We have also 


u{x)= +go 


1 _ llfi -^oll 

\x-Xo\\ r\\x-yi\\ 




_1 _ Il^~+1 J/all t = ®_ £ g _ r). 

\x-Xn\\ r\\x-yn+i\\J r 


Consequently, the conditions 2,3 satisfy too. 

5. The density of the equilibrium distribution 

The density of eqnilibrinm distribution on dB{xo,R) denote by 

Pr{x), xedB{xo,R) 

and the density of eqnilibrinm distribntion on dB{yo,r) denote by 

prix), xedB{yo,r) 

For Pr{x), X G dB{xo,R) we have 

1 /n ~ 11^0 — 


pR \x} = 




similarly for Pr{x), x G dB{yo,r) we have 

Pr\-^) A / . II ^ ^ ll'H 

47rr ^ |p-|/n|r 

Let us note that the potential fnnction for the eqnilibrinm distribntion 
eqnals 

U{x) = ^, xedB{xo,R) 

it 

U{x) = —, xedB{yo,r). 
r 
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If the point x is outside of the conductor E, i.e. x ^ E then we have: 
1 / PR(y] ,11 Pr 


U{x) = 


47ri?2 


'dB{xo,R) 


llx - y\\ 47rH F 


rds. 


6. The auxilary results 

Let Xn G B{xo,R) and Pm G B{yQ,r). Since 

^ Vm 


is harmonic function, so 

1 f R'^ - \\Xq - XnW^ 


r, X e B{xo, R) 


AtiR 


'dB{xo,R) 


\x — x„ 


x-yv 


- yn 


ds = 


Xn Vn 

\Xn I/m 


Now let us consider the case if bothe points are inside of the ball 
B{xq, R). Let Xn, Xm £ B{xo, R)- Remember that if ym-i is simmetric 
to Xm with respect the sphere dB{xo, R) then 


Since 


R\\x-Xm\\ = Po - ym-i\\\\x - ym-i\\, X E dB{xo,R). 


xeB{xo,R) 


Ill/m—1 x\ 

is harmonic function, so 

1 f R^ - \\Xo - Xn\\^ 


AnR 


'dB{xo,R) 


\X — Xr. 


R 


R^ — ||xo — X. 


X Xffi I 
2 


-ds = 

1 


11^0 I/m— l||47ri? \JdB{xo,R) ^n||^ ||^ |/m—1| 


ds I = 


R 


R 


\\Xn - ym-l\\\\Xo - ym-l\\ (Po “ “ Po “ Xr. 

_ ||:?0 - :?n||Po - Xm\\ 

R{R^ - ||fo - :?n||||:?o - Xm\\)' 
Applying the operator we obtain: 

1 f R^-\\xo-Xn\\^ 


I/m—11 


X Ytl X 


47ri? 


-ds = 


'dB{xo,R) 


X — Xr. 


Xm. - X 


Vo - Xo 


R\\Xn - foil 


||y 0 - foil - ||fo - fmllllfo - fnll)^' 
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7. The force of interaction 

The force acting on the sphere dB{xo, R) equals 

1 f 


Fr = - 


'dB{xo,R) 


WU{x)pR{x)ds = 


1 


Stt-R 

1 

Stt-R 


EE QnQr 

n=0 m=0 
oo oo 

'y ^ y ^ QmQn 


Xm-X - ||fo - Xn\\^ , , 
-^ 77 ^^- ds+ 


ldB{xo,R) 


X — Xr. 


\X — Xr. 


Vm-X R‘^ - \\Xo - XnW"^ _ 


m=0 n=0 


where 


^-^^EE 


'dB{xo,R) ~ ymir 

= p ^0 - Vo 
- foil’ 

QnQmR\\Xn ^o|| 


X - Xr, 


||XO-^m||||:?0-^n||)2 


EE 


QnQr. 


QnQm Ru{n) 


EE^ 


QnQr 


n=0 m=0 
oo oo 

^ ^ 2 (i ?2 _ u(m)u(n)y ^ ^ 2(d — uim) — v{n)Y' 

n= 0 m =0 ^ \ J \ JJ n= 0 m =0 ^ ^ ^ 

Similarly, The force acting on the ball B{yo,r) equals 


Fr = — 
2 


/ VU{x)pr{x)ds = F ? 

'dB{yo,r lll/Q —a^oH 


where 


^-EE 


rqnqmWyo - 


llfo-fn||||fo-fm||)2 


EE 


QnQr 


n=0 m=0 
oo oo 


EE^ 


Qnqn 


qnqmrv[m) 

^■^^2(d — u(n) — v(m)R ' ■^^2(r‘^ — v(n)v(m)y 

0—1) m=\\ ^ ^ ^ \ / / r)=n rr}=(\ ^ \ / \ / / 


EE 


n=0 m=0 ^ ^ ^ \ // ^_Q 

In numerical calculations it is useful to take into account the following 
inequality 

Jp QnQmRUyn) 

^ A—/ A—/ 9t — 11 (m^ii (n^R 


qnQr. 


^ “ u{m)u{n)y 2{d - u{m) - v{n)y 


OO OO 

s E E 


I Qn 11 Qn 


n=0 m=0 
oo oo 


< 


|Qn| |Qm|'^('^) 


< 


(d — u(n) — v(m)Y 2(r‘^ — u(n)u(m)Y 

n=N+lm=0 ' ' \ // ri=f)m=N+l ' \ / \ // 


< A 


Rr 


{d — R){d — r) 



So the approximation converges as geometric progression. Similar re¬ 
sults we have for Fr. 


8. Conclusion 

In some values of parameters, one of the ball vibrates. Note that in 
considered model the forces that keep the charges on the boundaries of 
the balls are not explicitly presented. So the system is not close and 
the Newton’s law can not satished. 
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9. Appendix. 


TwoBalls 

for n = 2:N-|-1; 

u(n) = r*r/(d - v(n-l)); 

v(n) = R*R/(d- u(n-l)); 

end 

Q1 = ones(l, N-l-1); 
ql = ones(l, N-|-l); 
for n = 2:N-|-1; 

Ql(n) = -ql(n-l)*u(n)/R; 
ql(n) = -Ql(n-l)*v(n)/r; 
end 

M = [sum(Ql(2:2:N)), sum(Ql(l:2:N)); 
sum(ql(2:2:N)), sum(ql(l:2:N));]; 



A = M/Q; 

QQ = A(l, 
qq = A(2, l).*ql; 

[U, V] = ndgrid(u, v); 

[ill, ml]=ndgrid(l:N+l,l:N+l); 

FRll = (QQ’*qq)./(2*(d - V - U).*(d - V - U)); 

FR12 = (QQ’*qq).*U.*R./(2*(R*R- U.*U’).*(R*R- U.*U’)); 
Fr21 = (qq’*QQ)./(2*(d - V - U).*(d - V - U)); 

Fr22 = (qq’*qq).*V.*r./(2*(r*r - V’.*V).*(r*r - V’.*V))); 

FR = sum(FRll(:))-sum(FR12(:)); 

Fr = sum(Fr21(:))-sum(Fr22(:)); 
testTwoBalls 
c =0.01:0.01:5; 

FR = Hi 

R' = D; 

for k = l:numel(c); 

[FR(k), Fr(k)] = TwoBalls(c(k)); 
end 

for k=l:numel(c); 

O(k)=0; 

end 

figure(l) k=30:l:nuniel(c); plot(k,FR(k),’r’,k,Fr(k),’b’,k,0(k),’y’) 
figure(2) k=l:l:30; plot(k,FR(k),’r’,k,0(k),’y’) 
figure(3) k=l:l:30; plot(k,Fr(k),’b’,k,0(k),’y’) 


8 




Figure 1. Let R=1 be the radius of the big ball and 
r=0.8 be radius of the small ball. Let d=R+r+ck be 
the distance of the ball’s centers, where c = 0.001, and 
A; = 1 : 17. In the picture you see Fr. 



Figure 2. Let R=1 be the radius of the big ball and 
r=0.8 be radius of the small ball. Let d=R+r+ck be 
the distance of the ball’s centers, where c = 0.001, and 
A; = 20 : 70. In the picture you see Fr. 
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Figure 3. Let R=1 be the radius of the big ball and 
r=0.8 be radius of the small ball. Let d=R+r+ck be 
the distance of the ball’s centers, where c = 0.001, and 
k = 500 : 2000. In the picture you see Fj. and Fr for big 
distences. 
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